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Abstract 

Background: Nasopharyngeal carcinoma (NPC) has a high metastatic feature. At/V'-Dinitrosopiperazine (DNP) is 
involved in NPC metastasis, but its mechanism is not clear. The aim of this study is to reveal the pathogenesis of 
DNP-involved metastasis. 6-1 OB cells with low metastasis are from NPC cell line SUNE-1, were used to investigate 
the mechanism of DNP-mediated NPC metastasis. 

Results: 6-1 OB cells were grown in DMEM containing 2 H 4 -L-lysine and l3 C 6 15 N 4 -L-arginine or conventional L-lysine 
and L-arginine, and identified the incorporation of amino acid by matrix-assisted laser desorption/ionization 
time-of-f light mass spectrometry. Labeled 6-1 OB cells were treated with DNP at 0 -18 |jM to establish the 
non-cytotoxic concentration (NCC) range. NCC was 0-10 |iM. Following treatment with DNP at this range, the 
motility and invasion of cells were detected in vitro, and DNP-mediated metastasis was confirmed in the nude 
mice. DNP increased 6-1 OB cell metastasis in vitro and vivo. DNP-induced protein expression was investigated using 
a quantitative proteomic. The SILAC-based approach quantified 2698 proteins, 371 of which showed significant 
change after DNP treatment (172 up-regulated and 199 down-regulated proteins). DNP induced the change in 
abundance of mitochondrial proteins, mediated the status of oxidative stress and the imbalance of redox state, 
increased cytoskeletal protein, cathepsin, anterior gradient-2, and clusterin expression. DNP also increased the 
expression of secretory AKR1B10, cathepsin B and clusterin 6-1 OB cells. Gene Ontology and Ingenuity Pathway 
analysis showed that DNP may regulate protein synthesis, cellular movement, lipid metabolism, molecular transport, 
cellular growth and proliferation signaling pathways. 

Conclusion: DNP may regulate cytoskeletal protein, cathepsin, anterior gradient-2, and clusterin expression, 
increase NPC cells motility and invasion, is involved NPC metastasis. 

Keywords: Dinitrosopiperazine, Carcinogen, Nasopharyngeal carcinoma, Metastasis, Quantitative proteomics 



Background 

Nasopharyngeal carcinoma (NPC) is a common malignant 
cancer in southern China [1]. Epidemiological investiga- 
tions have revealed that the incidence of NPC has 
remained high in endemic regions, particularly in south- 
ern China with an incidence of 30-80 per 100,000 people 
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per year [2]. NPC has the feature of high invasion and 
metastasis, and cervical lymphadenopathy is often the 
only clinical manifestation at initial diagnosis of NPC 
patients [3]. Therapeutic failure in advanced NPC has 
resulted from both high rates of local recurrence and dis- 
tant metastasis. 

In Chinese populations in high-incidence regions, the 
relative risk of NPC is related to their eating habits of the 
region, particularly with dietary intake of salt-preserved fish 
[2,4-6]. The process of salt preservation is inefficient and 
foods can become partially putrefied , consequently, these 
foods accumulate significant levels of nitrosamines [7,8], 
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which are known carcinogens [7,9,10]. A/jAf'-Dinitrosopi- 
perazine (DNP) is a predominant volatile nitrosamine in 
salted fish [11,12]. The carcinogenic potential of DNP in 
salt-preserved fish is supported by experiments in rats, 
which develop malignant nasal and NPC [13-15]. Further- 
more, DNP can induce malignant transformation of 
human embryonic nasopharyngeal epithelial cells [16]. Our 
previous works have shown that DNP induces rat NPC 
and shows organ specificity for nasopharyngeal epithelium, 
and found that DNP triggers over-expression of hot shock 
protein 70 and mucin 5B [17]. Additionally, DNP induces 
ezrin phosphorylation at Thr567 through activating Rho 
kinase and protein kinase C, and increases motility and in- 
vasion of NPC cells [18]. In the present study, to fully 
understand the mechanism of DNP-mediated NPC inva- 
sion and metastasis, we used a stable isotope labeling 
with amino acids in cell culture (SILAC) to further analyze 
the proteomic changes caused by DNP. We found that 371 
proteins were regulated by DNP, most of which were not 
previously reported to be involved in NPC metastasis. Ana- 
lysis of this vast information provides us with better under- 
standing of the complex regulatory mechanism of NPC 
high metastasis. Using bioinformatics analysis, we detected 
many novel signaling components in DNP-regulated sig- 
naling pathways. 

Methods 

Cell culture and stable isotope labeling 

NPC cell line 6-10B was derived from the cell line 
SUNE-1, and has a low metastatic ability [19]. Thus, 
6-10B cells were used in the present study to investigate 
DNP-mediated NPC metastasis. DNP was a carcinogens 
specially for NPC and its chemical structure is shown 
in Figure 2A. Additionally, heavy lysine and arginine 
( 2 H 4 -L-lysine and 13 C 6 15 N 4 -L-arginine) were purchased 
from Sigma- Aldrich. 6-10B cells were grown in DMEM 
containing 2 H 4 -L-lysine and 13 C 6 15 N 4 -L-arginine ("heavy") 
or conventional L-lysine and L-arginine ("light") supple- 
mented with dialyzed fetal bovine serum. After six cell 
doublings, we assessed the labeled amino acids in cells, 
and then identified whether cells were completely incor- 
porated by matrix-assisted laser desorption/ionization 
time-of-flight mass spectrometry using an Autoflex 
MALDI-TOF-MS instrument (Bruker Dalton). Mass 
spectra were searched against a database of human pro- 
teins and subsequently quantified using Mascot Server 
(www.matrixscience.com). 

3-(4,5-dimethylthiazol-2-yl)-5-(3- carboxyme -thoxyphenyl)- 

2- (4-ulfophenyl) -2H-tetrazolium assay 

To determine the non-cytotoxic concentration of DNP, 

3- (4,5-dimethylthiazol-2-yl)-5-(3- carboxyme -thoxyphe- 
nyl)- 2-(4-ulfophenyl) -2H-tetrazolium (MTT) assay was 
performed to determine the viability of "heavy" labeled 



6-10B cells. Briefly, "heavy" labeled 6-10B cells were seeded 
in 96-well plates at a density of 5 x 10 3 cells / well and trea- 
ted with DNP at a concentration between 0 and 18 \iM 
at 37°C for 24 h. Thereafter, 20 ul MTT (5 mg/ml, 0.5% 
MTT) was added per well for 4 h. The viable cell number 
per dish is directly proportional to formazan production, 
which can be measured spectrophotometrically at 492 nm 
following solubilization with isopropanol. 

Lactate dehydrogenase assay 

To further evaluate the non-cytotoxic concentration of 
DNP in "heavy" labeled 6-10B cells, lactate dehydrogenase 
(LDH) activity in cell culture media was detected after 
DNP treatment. Briefly, "heavy" labeled 6-10B cells were 
seeded in 6-well plates at a density of 2 x 10 4 cells/well 
and treated with DNP at a concentration between 0 and 
18 (iM at 37°C for 24 h. After the exposure period, media 
were collected for LDH activity measurement using the 
LDH assay kit (Autec Diagnostica). 

DNP treatment and protein preparation 

DNP crystals were dissolved in DMSO. Appropriate 
amounts of the DNP stock solution were added into the 
culture medium to achieve the indicated concentrations 
(DMSO concentration, 0.1%) and then incubated with 
cells for the indicated time periods. At approximately 
60% confluence, the "heavy" labeled 6-1 0B cells were 
treated with 10 uM DNP for 24 h according to MTT 
assay data, while the "light" labeled 6-10B cells were 
treated with only 0.1% DMSO, served as the control. 
The treated cells were then harvested and suspended 
with lysis buffer. Lysate was centrifuged at 13,200 rpm at 
4°C for 30 min. Supernatant fractions were collected and 
protein concentrations were determined using BCA 
assay kit (Pierce). 

Cell invasion and motility assay 

Cell invasion and motility were assayed according to 
methods described previously with minor modifications 
[18]. For the invasion assay, "heavy" labeled 6-10B cells 
were treated with the indicated concentrations of DNP 
for the indicated times. The treated cells were seeded 
into Boyden chamber with Matrigel (Neuro Probe, 
Cabin John, MD) at the upper part at a density of 1.5 
x 10 4 cells/well in 50 [A serum-free medium and incu- 
bated for 12 h at 37°C. The bottom chamber also con- 
tained standard medium with 20% fetal bovine serum. 
The cells invaded to the lower surface of chamber 
membrane were fixed with methanol and stained with 
hematoxylin and eosin. The invaded cell numbers were 
counted under a light microscope. The motility assay 
was performed as described in the invasion assay with- 
out Matrigel coating. 
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Figure 1 Non-cytotoxic concentration of DNP in stable isotope- 
labeled 6-1 OB cells. A, Structure of DNP, an N-nitroso compound. 
B, Stable isotope-labeled 6-1 OB cells were treated with the indicated 
DNP concentration and then subjected to IVuT cell viability assay as 
described in the Materials and Methods. "OD" indicates the relative 
optical density at 492 nm. C, After stable isotope-labeled 6-1 OB cells 
were treated with the indicated DNP concentration, and cell media 
were subjected to the LDH assay. LDH activity is per 1 L of media. 
Data are presented as the means ± standard deviation from three 
independent experiments, statistically analyzed using Student's 
t- test (*, p < 0.05; a=0.05). 



Evaluation of the effect of DNP on NPC metastasis 
in nude mice 

Nude mice experiments were performed as previously 
described [18]. Twenty BABL/c nude mice (approxi- 
mately 5-6 weeks old) were purchased from the Animal 
Center of Central South University. All animal studies 
were conducted according to the standards established 
by the Guidelines for the Care and Use of Laboratory by 
Animals of Central South University. Additionally, the 
present study protocols were approved by the ethical 
committee at Central South University. Briefly, 100- ul 
aliquots of 6-10B cell suspensions (1 x 10 4 cells) were 
mixed with Matrigel and injected respectively into the 
tail veins of the 20 nude mice. They were then randomly 
divided into two groups, DNP-treated and control 
groups, containing 10 mice per group. The DNP-treated 
group was abdominally injected with DNP at a dose of 
40 mg/kg (body weight) twice a week for 60 days using a 
1-ml sterile syringe. The control group was treated with 
0.1% DMSO. After DNP treatment, the metastasis of 
6-10B cells to the lung, liver, and lymph nodes was 
observed. Their metastatic abilities were evaluated by 
counting tumor metastatic foci on day 60 after the 
injection. 

Gel electrophoresis and in-gel trypsin digestion 

Prior to gel electrophoresis, equal amounts of DNP- 
treated and untreated cell proteins were mixed, separated 
using 10% SDS-PAGE (4 - 12% Bis-Tris Novex minigel, 
Invitrogen), and stained silver solution to visualize the 
gel bands. The entire protein gel lanes were horizontally 
excised and cut into 48 slices each, and then destained, 
reduced, alkylated and digested with gold-trypsin at 
37°C overnight as described previously [20]. The result- 
ing tryptic peptides were extracted by 90% acetonitrile 
(Fisher) and 2.5% trifluoroacetic acid (Promega), lyophi- 
lized in a SpeedVac, and dissolved in 1% formic acid and 
2% acetonitrile before liquid chromatography-tandem 
mass spectrometry (LC-MS/MS) analysis. 

LC-MS/MS analysis 

The peptide mixtures were separated using Finnigan 
Surveyor high-performance liquid chromatography system 
(Thermo Electron, San Jose, CA) on a C18 reverse phase 
column, which was coupled online to a linear ion trap/ 
Orbitrap (LTQ-Orbitrap) mass spectrometer (Thermo 
Electron, San Jose, CA). Briefly, the peptide mixtures 
were first loaded onto a C18-reversed phase column 
(100-um inner diameter, 10-cm long, 3-um resin from 
Michrom Bioresources, Auburn, CA), and then sepa- 
rated at a maximal flow rate of 300 nl/min controlled by 
IntelliFlow technology. The peptide mixtures were sepa- 
rated using the following parameters: 1) mobile phase A: 
0.1% formic acid, 2% the acetonitrile, Dissolved in water; 



Li et al. BMC Biochemistry 2012, 13:25 
http://www.biomedcentral.eom/1 471-2091/1 3/25 



Page 4 of 13 



2) mobile phase B: 0.1% formic acid, dissolved in aceto- 
nitrile; 3) flow rate: 300nl/min; 4) gradient: B-phase 
increased from 5% to 35%, 120min. Next, the eluate was 
online analyzed online in LTQ-Orbitrap mass spec- 
trometer operated in a data-dependent mode, the 
temperature of the heated capillary was set to 200°C, 
and the spray voltage was set to 1.85 kV. Full-scan 
MS survey spectra (m/z 400-2,000) in the profile 
mode were acquired in Orbitrap with a resolution of 
60,000 at m/z 400 after the accumulation of 1,000,000 
ions, and followed by five MS/MS scans in LTQ with 
the following Dynamic Exclusion settings: a repeat 
count of 2, a repeat duration of 30 s, and an exclusion 
duration of 90 s. The lock mass option was enabled 
for survey scans to improve mass accuracy [21]. The 
data were acquired using Xcalibur (Thermo Electron, 
version 2.0.7). 

Protein identification, quantification and bioinformatics 
analysis 

Protein identification and quantification were performed 
as previously described with minor modifications [22,23] . 
Briefly, the mass spectrometric raw data were analyzed 
using MaxQuant 1.0.13.13 software and the derived peak 
lists were searched using the Mascot search engine 
(Matrix Science, version 2.2.04, London, UK) against a 
concatenated real and false International Protein Index 
human protein database (V3.52). Mascot search results 
were further processed by MaxQuant 1.0.13.13 at the 
false discovery rate of 1% at both the protein, peptide, 
and site levels. The normalized heavy versus light (H/L) 
ratios, significance, and variability (%) were automatically 
produced by MaxQuant 1.0.13.13 software. The final 
reported protein ratio represents a normalized ratio of 
H/L SILAC obtained in all technological repeats where 
the same protein was identified. International Protein 
Index numbers of all significantly regulated proteins 
and some unaltered proteins were imported into the In- 
genuity Pathway Analysis software tool (www.ingenuity. 
com) for bioinformatics analysis based on published 
reports and databases such as Gene Ontology, Uniport, 
and TrEMBL. 

Western blotting analysis 

Western blotting was used to validate the expression 
levels of eight dysregulated proteins in DNP-treated and 
untreated 6-10B cells as described above. 6-10B cells 
were treated with 5, 10, 20 uM for dose-course and trea- 
ted with 10 uM for 6, 12, 18, 24, 36, 48 h for time- 
course. After treatment, supernatants were centrifuged 
at 300 x g for 4 min and 2000 x g for 8 min to remove 
dead cells and cell fragments, and proteins were concen- 
trated by centrifugal ultrafiltration using Microcon YM-3 
Centrifugal filters (Millipore, Billerica, MA, USA). The 



treated cells were disrupted with 0.6 ml lysis buffer 
[1 x PBS, 1% Nonidet P-40, 0.1% SDS, and freshly added 
100 (ig/ml PMSF, 10 (ig/ml aprotinin, 1 mM sodium 
ortho vanadate]. Cell lysates were then subjected to cen- 
trifugation of 10000 x g for 10 min at 4°C. Equal protein 
amounts of cell lysates and culture supernatants were 
separated by 10% polyacrylamide gel electrophoresis, 
and transferred onto nitrocellulose membranes (Bio-rad). 
The membranes were subsequently incubated with 5% 
non-fat milk in Tris-buffered saline containing 0.05% 
Tween-20 for 1 h to block non-specific binding and then 
overnight with antibodies against aldo-keto reductase 
(AKR) 1B10, S100P, cathepsin B, cathepsin D, ferritin, 
a-E-catenin (Cell Signaling Technologies), or clusterin, 
AGR2, and GAPDH (Santa Cruz.), then incubated with 
the secondary antibody for 1 h at room temperature. The 
band signal was developed using 4-chloro -1-napthol/ 
3,3-o-diaminobenzidine, and relative photographic density 
was quantitated using a gel documentation and analysis 
system (Pierce, Rockford, USA). 

Gene transfect and wound-healing assays 

Wound-healing assay was performed as previously 
described with minor modifications [24]. 6-10B cells 
(2 x 10 6 ) were seeded in 10-mm plates at 37°C for 24 h, 
and transiently transfected with si-AGR2 or si-mock 
(Dharmacon) [25] using Lipofectamine 2000 reagent (Life 
Technologies, Inc.) following the manufacturer s suggested 
protocol, and then confluent monolayer of the transfected 
cell was wounded using a plastic tip. Cells were treated 
with DNP at 10 uM, and then photographed after 48 h. 
The cells moving cross the boundaries lines were counted. 
The transfect cell samples were harvested, and total 
proteins were extracted. These protein samples were sub- 
jected to Western blotting analysis. 

Results and discussion 

In this study, quantitative proteomics with SILAC were 
used to identify the different protein of 6-10B cells with or 
without DNP treatment. As the first step 6-10B cells were 
labeled with amino acid, and then we assessed the incorp- 
oration efficiency of 2 H 4 -L-lysine and 13 C 6 15 N 4 -L-arginine 
in 6-10B cells for full incorporation in all proteins after six 
cell doublings. Three peptides, VEVTEFEDIK (Figure 1A), 
GHYTEGAELVDSVLDWR (Figure IB) and LRQPFFQK 
(Figure 1C) were separated by 4 Da, 10 Da, and 14 Da cor- 
responding to the mass difference between the above light 
and heavy isotopes. The entire signal corresponded to the 
heavy peptide, indicating that incorporation of 2 H 4 -L-lysine 
or 13 C 6 15 N 4 -L-arginine was complete. To illustrate the 
quality of the protein identifications reported, we present 
MS and MS/MS spectra of clusterin and AKR1B10 from 
the data obtained from the LTQ-Orbitrap mass spectrom- 
eter (Figure ID, E). 
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Figure 2 Labeling of 6-1 0B cells with 2 H 4 -L-lysine and 13 C 6 15 N 4 -L-arginine. 6-1 0B cells were cultured in DMEM containing 2 H 4 -L-lysine and 
13 C 6 15 N 4 -L-arginine or conventional L-lysine and L-arginine as described in the Materials and Methods. A, the peptide VE\^EFEDIK has a charge 
state of 1 + and contains one lysine; B, the peptide GHYTEGAELVDSVLDWR has a charge state of 1 + and contains one arginine. C, the peptide 
LRQPFFQK has a charge state of 1 + and contains one lysine and one arginine. The peaks marked with blue and red represent the light isotope 
signal and the heavy isotope signal, respectively. Representative MS/MS and MS spectra from clusterin (D) and AKR1B10 (E). The inset shows the 
relative ratio of heavy to light versions of each precursor ion. 



DNP is a carcinogenic A/-nitroso compound, and its 
chemical structure is shown in Figure 2A. Although it is 
known that the non-cytotoxic concentration of DNP to 
6-10B cells is 0-4 |iM, even up to 6 [iM [18], we deter- 
mined the non-cytotoxic concentration of DNP in stable 
isotope-labeled 6-10B cells, cell viability was not signifi- 
cantly altered at 0-10 (iM DNP compared with control 
cells (0.1% DMSO; Figure 2B, *, p < 0.05). To further 
confirm that 0-10 \iM DNP was non-cyto toxic, LDH ac- 
tivity in the cell culture media was detected after DNP 
treatment. The data revealed that LDH activity was not 
significantly altered by treatment with DNP concentra- 
tions between 0 and 10 \iM (Figure 2C, *, p < 0.05). Thus, 
in all subsequent experiments, the concentration of 
10 \iM DNP was used. 

Although previous work has shown that DNP is 
involved in NPC metastasis, we first confirmed here 
that DNP mediated NPC metastasis. A Matrigel- coated 



Boyden chamber was used to measure invasion. 6-10B 
cells were treated with DNP at 0, 2, 4, 6, 8, and 10 [iM 
for 24 h and then seeded into the Boyden chamber. 
The cells that invaded the lower chamber were counted. 
The invaded cells increased dose-dependently after DNP 
treatment (Figure 3A-c, lanes 4 to 6 vs lane 1; *, p < 0.05). 
Compared with the control, the increase was 4.1 -fold with 
8 \iM DNP (Figure 3A-c, lane 5). For detecting 6-10B cell 
motility with DNP treatment, the treated cells were seeded 
into a Boyden Chamber uncoated with Matrigel, and 
motile cells were counted. A similar effect was observed 
for the motility of DNP-treated cells (Figure 3B-c, lanes 4 
to 6 vs lane 1; *, p < 0.05). The cell motility increased by 
5.6-fold after treatment with 8 [iM DNP (Figure 3B-c, 
lane 5). To further confirm DNP -involved metastasis 
in vivo, the treated 6-10B cells were mixed with Matrigel, 
and then were injected into the tail veins of BABL/c 
mice. Tumor metastatic nodes of 6-10B cells in the 
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Figure 3 DNP-mediated 6-1 0B cell invasion and motility in vitro or metastasis in vivo. A, stable isotope-labeled 6-1 OB cells were treated 
with the indicated DNP concentration for 24 h. The treated cells were subjected to analyses for motility and invasion as described in the Materials 
and Methods. A, invasion of 6-1 OB cells at various concentrations. Arrow, invaded cell. B, motility of 6-1 OB cells at various concentrations; Arrow, 
motile cell. The data were statistically analyzed by one-way analysis of variance with post hoc Dunnett's test (*, p < 0.05). Scale bar, 5 uJVL 
C, twenty nude mice were injected with 6-1 OB cells in Matrigel through the tail vein ( 1 x 1 0 4 cells/mouse), and then randomly divided into two 
groups with 10 mice per group. The DNP group was abdominally injected with DNP using a stumped needle at a dosage of 40 mg/kg 
(bodyweight), twice a week for 60 days. The control was injected with 0.1% DMSO. After 60 days, tumor metastatic foci of 6-1 0B cells were 
observed. Scale bar, 10 uM. Arrow, tumor metastatic foci. 



lungs, livers and lymph nodes were detected. Metastatic 
foci in mice lungs were significantly observed in nude the 
mice with DNP treatment (Figure 3C, left panel vs right 
upper line), and pathologically confirmed under micro- 
scope (Figure 3C, left panel vs right down line). These 
data indicated that DNP mediates NPC metastasis in 
vitro and in vivo. 

To fully investigate the mechanism of DNP-mediated 
NPC metastasis, SILAC coupled with LC-MS/MS was 
used to identify and quantify the proteomic differences. 
A total 2853 proteins were detected, and 2698 (94.57%) 
proteins could be quantified. Of these 2698 protein, 172 
were calculated to be highly up-regulated, and 199 were 
significantly down-regulated at a ratio H/L >2.0 or ratio 



H/L <0.5 and p < 0.05 (Figure 4A). To gain functional 
insight into the cellular proteome, the 172 up-regulated 
and 199 down-regulated proteins were respectively 
assigned to different molecular functional classes and 
subcellular annotations according to the underlying bio- 
logical evidence from the Gene Ontology database 
(false discovery rate < 0.05). Because some proteins 
generally have more than one component annotation 
or function annotation, the sum of each category may 
be higher than 100%. The 15 most abundant terms 
are shown in Figure 4B, C, with additional data shown 
in Additional file 1: Table SI. Mitochondrion proteins 
and proteins related to junctional mechanisms were 
highlighted in up-regulated and down-regulated 
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Figure 4 Proteome quantitation, significance, and classification analysis. A, signal intensities of all quantified proteins after DNP treatment 
are shown as a function of their fold change. Ratios of most proteins distributed around 1.0 (Log 10 (ratio) = 0), indicating that whole proteins of 
the two groups of cells were mixed equally and 6-1 OB cells were fully labeled. The spread of the cloud was lower at high abundance, indicating 
that quantification is more precise, and their fold change level is indicated in blue, red, and green, respectively. To reduce test error, proteins at 
more than 2.0-fold or less than 0.5-fold and p < 0.05 were deemed to indicate significantly changed proteins induced by DNP treatment. 
Together, this gives quantifiable results for 371 proteins being significantly altered upon DNP, with 172 and 199 proteins being up- and down- 
regulated, respectively. The 15 most abundant terms of component (B) and function annotation in disregulated proteins (C) are shown. As each 
protein is generally assigned to more than one term, the percentage of proteins in each term is shown instead of the total number to avoid 
redundancy. Number distribution of the 15 most abundant terms for up-regulated or down-regulated proteins (D). Oxidoreductase activity-1, 2, 3, 
and 4 represent respectively oxidoreductase activity, oxidoreductase activity acting on NADH or NADPH, oxidoreductase activity acting on the 
CH-CH group of donors, and oxidoreductase activity acting on the CH-CH group of donors, NAD or NADP as acceptor (E). 



proteins individually (Figure 4D), suggesting that fur- 
ther exploration at subcellular levels is necessary. 
Functional analysis of these differential proteins 
showed that DNP-treated high metastatic 6-10B cells 
demonstrated significant changes in oxidoreductase 
activity, cofactor binding, and cytoskeletal protein 
binding (Figure 4E). 

Proteins that changed significantly in DNP-treated cells 
were mapped to 15 specific functional networks with 
each network containing 11 or more "focus" members 
(Figure 5 A, Additional file 2: Table S2). The four net- 
works of interest correspond to the following: (A) Cancer, 



Renal and Urological Disease, Cell Death (Figure 5B); 
(B) Cancer, Reproductive System Disease, Cell Death 
(Figure 5C); (C) Cellular Movement, Lipid Metabolism, 
Molecular Transport (Figure 5D), and (D) Protein Syn- 
thesis, Cell Death, Cellular Growth and Proliferation 
(Figure 5E). Proteins that are present in these pathways 
and that were identified in our analysis as up-regulated are 
depicted in red, and proteins that were identified as down- 
regulated are shown in green. Proteins known to be in 
the network but that were not identified in our study are 
depicted in white. The shade of the color indicates the 
magnitude of the change in protein expression level. 
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(See figure on previous page.) 

Figure 5 Ingenuity Pathway Analysis of proteins induced by DNP. A, overview of 15 specific functional networks, each of which contains 1 1 
or more "focus" proteins (proteins that were significantly up- or down-regulated). Each box contains an arbitrary network number. B, cancer, renal 
and urological disease, and cell death. C, cancer, reproductive system disease, and cell death. D, cellular movement, lipid metabolism, and 
molecular transport. E, protein synthesis, cell death, and cellular growth and proliferation. Red, up-regulated proteins; Green, significantly down- 
regulated proteins; White, proteins known to be in the network but were not identified in our study. The color depth indicates the magnitude of 
the change in the protein expression level. Lines connecting the molecules indicate molecular relationships. Dashed lines indicate indirect 
interactions, and solid lines indicate direct interactions. The arrow styles indicate specific molecular relationships and directionality of the 
interaction. Abbreviations are shown in Additional file 2: Table S2. 



To confirm the SILAC results, eight proteins with dif- 
ferent fold changes, AKR1B10, clusterin, cathepsin B, 
cathepsin D, ferritin, a-E-catenin, AGR2, and S100P 
were chosen to validate SILAC results. Western blotting 
results showed that the ratios of eight representative 
proteins between treated and untreated cells showed 
either close-degrees or similar fold changes consistent 
with those obtained from SILAC (Figure 6A). Quantifica- 
tion results by SILAC of the eight proteins are shown in 
Additional file 3: Table S3. Some of them are secretory 
proteins, to determine whether DNP also induces these 
secretory proteins, AKR1B10, cathepsin B and clusterin 
in the culture supernatants of DNP treated 6-10B cells 
were detected. The results showed that AKR1B10, 
cathepsin B and clusterin dramatically increased in the 
culture supernatants after DNP treatment (Figure 6B). 
These findings imply that DNP may induce 6-1 OB cells 
to secrete AKR1B10, cathepsin B and clusterin. 

To further confirm whether the different proteins involve 
DNP-mediated NPC metastasis, we chose high-expressed 
protein AGR2 as target. DNP induced AGR2 expression at 
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Figure 6 Differential protein expression confirmed using 
Western blotting. A, eight representative proteins, aldo-keto 
reductase (AKR) 1B10, clusterin, cathepsin B, cathepsin D, ferritin, 
a-E-catenin, AGR2, and S100P were detected in 6-1 OB cells with or 
without DNP treatment using Western blotting. GAPDH served as a 
loading control. The relative photographic density was quantitated 
using a gel documentation and analysis system. B, AKR1B10, 
Cathepsin B and Clusterin were detected in the cell culture 
supernatant of 6-1 OB cells with or without DNP treatment using 
Western blotting. WB, Western blotting. SILAC, stable isotope 
labeling with amino acids in cell culture. 



dose- and time- course (Figure 7A). As a metastasis- 
associated protein, AGR2 may play an important role in 
DNP-mediated metastasis. The next step is to observe 
DNP-mediated metastasis when AGR2 blocked. We used 
si-AGR2 to knockdown AGR2 (Figure 7B), and then used 
wound-healing assays to detect the cell motility of 
6-10B-siAGR2 with DNP treatment. Following si-AGR2 
transfect, DNP-mediated motility decreased, and conse- 
quently the cells were unable to migrate into the wound 
(Figure 7C, panel b vs. d and Figure 7D, lane 2 vs. 4). 
Hence, we concluded that AGR2 plays an important role 
in DNP-mediated metastasis. 

Conclusion 

In clinic, NPC has the features of high invasion and 
metastasis, but its mechanism has been unclear. As 
one of three carcinogen factors for NPC, the Epstein- 
Barr virus (EBV) has been proven to be involved in NPC 
metastasis through latent membrane protein 2A inducing 
epithelial-mesenchymal transition (EMT), however latent 
membrane protein is positive at only a 56.7% rate [26]. 
Recently, another important carcinogen factor, DNP was 
also found to be involved in NPC metastasis [18]. In the 
present study, using SILAC and a systematic data ana- 
lysis method, we obtained unbiased interpretation of 
NPC cell metastasis induced by DNP. Approximately 
2698 proteins were quantified and 371 of these proteins 
showed apparent alterations in expression levels after 
DNP treatment, involving the regulation of biosynthesis 
and energy metabolism, as well as cell adhesion or inva- 
sion. We speculated that biosynthesis, energy metabolism 
and invasion are associated with NPC metastasis 
mediated by DNP. Based on subcellular and biological 
function analysis, many differential proteins in the 
present study were located in mitochondrion, such as 
mitochondrial membrane part, and mitochondrial re- 
spiratory chain. Additionally, tumor cells with mitochon- 
dria damage or dysfunction were reported to enhance 
anti-apoptosis ability and invasion [27,28]. This suggests 
that mitochondrial dysfunction may be linked to metas- 
tasis of DNP-treated 6-10B cells. 

In the differential proteins mediated by DNP, oxidore- 
ductase activity and oxidoreductase activity acting on 
NADH or NADPH, the CH-CH group of donors, and 
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Figure 7 DNP-mediated 6-1 OB cell motility through AGR2. A, 6-1 OB cells were treated with the indicated concentration DNP for dose-course, 
and treated with 10 uM DNP for the indicated time for time-course. AGR2 expression in the DNP-treated cells was detected using Western- 
blotting. 6-1 OB cells (2 x 10 6 ) were transiently transfected with si-AGR2 or si-mock, incubated at 37°C for 16-24 h. Confluent monolayer of the 
transfect cells was then wounded using a plastic tip, and then treated with 10 uM DNP. (a) 6-1 OB cells treated with 0.1% DMSO. (b) 6-1 OB cells 
treated with 10 uM DNP. (c) The transfect with AGR2 si-RNA treated with DMSO. (d) The transfect treated with DNP. B, AGR2 expression in the 
above transfect cells with DMSO or DNP treatment using western-blotting. C, The treated cells were observed under microscope, and the cells 
migrating cross the boundaries lines in the center of the wells were counted. D, Numbers of cells that moved cross the lines (10 fields). The data 
are represented as the mean ± SD from three independent experiments. Results were statistically analyzed using one-way analysis of variance 
(ANOVA) with a post hoc Dunnett's test (* p < 0.05). The error bars represent SDs. Scale bar, 10 uM. Arrow, cells moving cross the boundary. 



the CH-CH group of donors, NAD or NADP as the 
acceptor related to proteins accounted for a large pro- 
portion. Peroxiredoxins 3, NADH-dehydrogenase ubi- 
quinone iron-sulfur protein 3 (NDUFS3), NADH- 
dehydrogenase ubiquinone 1 beta subcomplex subunit 
8 (NDUFB8), pirin, ferritin heavy chain, and AKR1 were 
significantly up-regulated in the high metastatic 6-1 0B 



cells with DNP treatment. Oxidative stress have been 
shown to play important roles in tumorigenesis and 
progression of tumors [29], in which there is aberrant 
or improper regulation of the redox status. The balance 
of redox state affects many physiological and patho- 
physiological processes of cells, its mechanisms include 
gene transcription, cell signal transduction, activity of 
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enzymes and biological macromolecules, cell prolifera- 
tion, adhesion, and apoptosis. These findings suggest 
that the significant change of oxidoreductase activity in 
high metastatic 6-10B cells with DNP treatment is cor- 
related with the status of oxidative stress and imbalance 
of the redox state. 

Cytoskeleton has been identified as a major target for 
destruction during apoptosis and is important under 
pathological conditions such as cancers [30]. The differ- 
ential proteins were distributed in the cytoskeleton, in- 
cluding N-myc downstream-regulated gene 1 protein, 
paxillin, and syntenin-1. Conversely, some proteins asso- 
ciated with the cytoskeleton were up-regulated, such as 
catenin alpha- 1, radixin, macrophage-capping protein, 
integrin beta-5, tubulin-specific chaperone D, tubulin 
beta 2C (TUBB2C), tubulin beta 2A, and tubulin 5 beta. 
And subcellular localization of these differential proteins 
is related to junctional mechanisms. Based on these data, 
we speculate that in high metastatic 6-10B cells with 
DNP treatment, dynamic modifications and remodeling 
in the cytoskeleton exist, and the dynamic alteration 
affects endocytosis, cell shape, cell motility, cell adhesion 
and invasion. 

Additionally, some important proteins directly related 
to metastasis were discovered in our study, such as, 
annexin A6, S100P, S100A4, hot shock protein 90B1, 
ferritin heavy chain, TUBB2A, and anterior gradient-2 
(AGR2, Additional file 3: Table S3). Cathepsin B, 
AKR1B10 and custerin were not only up-regulated in 
6-10B cells with DNP treatment, but also in the cell 
culture supernatant. Cathepsins, initially described as 
intracellular peptide hydrolases, play a role in invasion 
and metastasis of cancer [31]. In the present study, 
cathepsins B and D were respectively up-regulated 
7.9-fold (Additional file 3: Table S3) and 4.6-fold 
(Additional file 3: Table S3), respectively. Cathepsin B is a 
key enzyme in invasion and metastasis of malignant 
tumors. It is up-regulated in laryngeal cancer [32], cer- 
vical cancer [33,34], and bladder cancer [35], and its ex- 
pression level is correlated with metastatic potential. 
Cathepsin D, a lysosomal aspartate proteolytic enzyme 
that is similar to cathepsin B, also plays an important 
role in invasion and metastasis of cancer. It is up- 
regulated in metastasis of some malignant tumors, in- 
cluding primary laryngeal cancers correlated with neck 
lymph node involvement [36], gastric cancer with lymph- 
atic and/or blood vessel invasion [37], and breast cancer. 
Furthermore, Cheng et al. [38] found that significant 
cathepsin D expression occurred in lymph node metasta- 
sis versus primary NPC and was significantly correlated 
with advanced clinical stage, recurrence, and lymph node 
and distant metastasis. AGR2 was reported to be linked 
with several human cancers and induced metastasis 
[39]. Additionally, Dumartin, et al. [25] found that 



cathepsins B and D are downstream functional mole- 
cules of the proinvasive AGR2 in vitro, and AGR2, 
cathepsins B and D were considered to be essential for 
dissemination of pancreatic cancer cells in vivo. High 
expressed-cathepsins B and D in DNP-treated 6-10B may 
be mediated by AGR2, but it is also possible that DNP dir- 
ectly mediated cathepsins B and D. Additionally, DNP- 
induced 6-10B motility decreased when AGR2 blocked 
(Figure 7). We speculated that cathepsins B, D and AGR2 
expression mediated by DNP and AGR2 regulating cathe- 
psins B, D are involved in NPC metastasis. 

Significantly, AKR1 proteins were predominantly 
up-regulated in high metastatic DNP-treated 6-10B cells, 
including AKR1C1, AKR1B10, AKR1C3, and AKR1B1 
(Additional file 3: Table S3). Family members of AKR1C 
play a pivotal role in maintaining steroid homeostasis 
and catalyzing reductive detoxification of reactive alde- 
hydes and ketones, which are produced as a result of 
oxidative stress [40,41]. AKR1B10 is also correlated posi- 
tively with tumor size and lymph node metastasis [42], 
These findings suggest that DNP would affect oxidative 
stress and steroid homeostasis in 6-10B cells through 
the above aldo-keto reductase family 1 proteins, thereby 
increasing 6-10B cell metastasis. 

Higher clusterin levels were expressed in various malig- 
nant tumors with metastasis including ovarian [43], breast 
[44], and gastric cancers [45]. An emerging query, cluster- 
ins enhanced cell invasion and metastasis of tumors 
through EMT. Lee, et al. [46] found that clusterin was 
involved in Smad2/3 stability at the protein level, and 
believed that clusterin regulates transforming growth 
factor-beta signaling pathway by modulating the stability of 
Smad2/3 proteins and mediates EMT. Lenferink et al. [47] 
also found that clusterin gene expression was highly up- 
regulated throughout transforming growth factor-beta, and 
speculated that secreted clusterin served as an important 
extracellular promoter of EMT. In the present study, pro- 
teins related to EMT and cell adhesion were also dysregu- 
lated, including clusterin myosin- VI, catenin alpha- 1 
(CTNNA1), fibronectin type III domain-containing protein 
3B (FNDC3B, 2.1-fold), LI cell adhesion molecule 
(LI CAM), desmoplakin, plakophilin-3 (Additional file 3: 
Table 3), implying that the mobility of DNP-induced 6-10B 
cells is probably related to EMT and cell adhesion. 

The SILAC technique was used to conduct a comparison 
of the proteomes of 6-10B cell metastasis induced by DNP. 
A cooperative response, including many proteins, and a 
group of pathways were identified and some interesting 
clues were provided. DNP may induce a change in abun- 
dance of mitochondrial proteins, mediate the status of oxi- 
dative stress and the imbalance of the redox state, and 
increase cytoskeletal protein, cathepsin, AGR2, and clus- 
terin expression, and finally promote cell metastasis. DNP 
may be involved in NPC metastasis through regulation of 
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cancer protein synthesis, cellular movement, lipid metabol- 
ism, molecular transport, cell death, and cellular growth 
and proliferation signaling pathways. DNP may also induce 
6-10B cells to secrete AKR1B10, cathepsin B and clusterin. 
These dataset provide important clues for investigation on 
high metastatic NPG 

Additional files 



Additional file 1: Table SI. Gene Ontology analysis. 

Additional file 2: Table S2. Ingenuity Pathways analysis. 

Additional file 3: Table S3. Information for up- and down- regulated 
proteins identified in DNP-induced cell. 
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